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Quasiperiodic oscillations with frequencies in the range 5-50 Hz have recently 
been discovered in X-rays from two bright galactic-bulge sources and Sco X-l. 
We propose that these sources are weakly magnetic neutron stars accreting from 
disks in which the plas&n is clumped. The interaction of the magnetosphere 
with clumps in the inner disk causes oscillations in the X-ray flux with many 
of the properties observed. 


Quasiperiodic oscillations have recently been discovered in X-rays from 
the bright galactic sources GX 5-1, 1> 2 , Cyg X-2, 3. 4 anc | geo X-l. 5,6 The 
frequency fp of the oscillations varies from 20 to 36 Hz in GX 5-1, 28 to <13 
Hz in Cyg X-2 , and from 6-17 Hz in Sco X-l. In all the available data on GX 
5-1 and Cyg X-2, fp shows a strong positive correlation with the 1-10 keV 
count rate. The FWHM Afp of the peak in averaged power spectra varies from 4 
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to 13 Hz in GX 5-1 and from 12 to 20 Hz in Cyg X-2, depending on the count 
rate, with relative widths Af p /f p in the range 0.21-0.60. For Sco X-l, only 
the 5-15 keV count rate is available; f p is strongly positively correlated 
with this count rate during 10-20 minute quiescent intervals between flaring 
episodes but appears to vary erratically during transitions to extended low 
states. The FWHM of the peak in the Sco X-l power spectrum varies from 1 to 
10 Hz. All three sources show substantial red noise. In GX 5-1 the rms 
fractional intensity variation caused by the quasiper iodic oscillations and 
that caused by the red noise both remain approximately constant as the source 
Intensity increases, except perhaps at the highest intensities observed where 
the variation may be less. In Cyg X-2, on the other hand, the fractional 
variation caused by the oscillations remains approximately constant but that 
caused by the red noise increases as the intensity increases. Finally, in Sco 
X-l the red noise down to the lowest frequencies so far explored is relatively 
weak when the oscillations are clearly present but becomes stronger when the 
source is most intense. 

Here we suggest that these sources are weakly magnetic neutron stars 
accreting from Keplerian disks in which the plasma is clumped. Interaction of 
the magnetosphere with clumps in the inner disk causes the mass accretion rate 
and hence the X-ray Intensity to vary at a harmonic of the beat frequency 
given by the difference between the rotation frequency of the star and the 
rotation frequency of the inner disk. [This idea has been suggested 
previously as a possible explanation for some of the quasiper iodic 
oscillations observed in cataclysmic variables (ref. 7; D . Lamb and Warner, 
private communication)]. The low harmonics of this beat frequency are 
comparable to the oscillation frequencies observed in GX 5-1, 8,9 Cyg X-2, and 
Sco X-l for neutron star spin rates - 100 Hz, dipole magnetic fields - 10 9 G, 
and accretion rates ~ 10 18 g s - *. Such spin rates and dipole magnetic field 
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strengths are consistent with the hypothesis that systems like these are 
progenitors of the millisecond rotation-powered pulsars. 8 In addition, for 
neutron stars rotating near their equilibrium spin rates the low harmonics of 
t.iis beat frequency vary with accretion rate in a manner similar to the 
observed variations of the oscillation frequency with X-ray intensity. 9 

We show that the power spectra given by a simple version of this model 


are very similar to 

the 

spectra reported in 

GX 5-1 

and 

Cyg X-2. In 

particular, red noise 

at 

frequencies below 

that of 

the 

quasiperiodic 


oscillations is a logical consequence of the model. If all clumps in the 
inner disk contribute to the oscillations and if the properties of the clumps 
and the nature of their interaction with the magnetosphere remain unchanged, 
the strength of the red noise is proportional to that of the quasiperiodic 
oscillations, as observed in GX 5-1. If not, the strength of the red noise 
may increase even as the strength of the oscillations decreases, as observed 
in Cyg X-2 and Sco X-l. We present calculations which show that X-ray 
pulsations caused by beaming and rotation of the neutron star are strongly 
suppressed relative to quasiperiodic oscillations produced according to the 
present model, under the conditions thought to exist in bright galactic bulge 
sources. Several observational tests of the model are described. 

Physical Model 

The plasma flow pattern around a weakly magnetic neutron star accreting 
from a Keplerian disk at near the Eddington critical rate is expected to be 
similar to that around supermassive black holes in galactic nuclei accreting 
at near -critical races. 10-12 The inner disk is radiation-pressure-dominated, 
unlike the inner disks surrounding canonical accretion-powered pulsars with 
dipole magnetic fields -10 12 G. The luminosity of the inner disk is supplied 
in part by release of gravitational binding energy as the disk plasma spirals 
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In, and in port by radiation from the neutron star. The pressure of the 
radiation leaving the inner disk alone implies a disk height h > 1.5(M/Mg)R, 
where Mg is the mass flux corresponding to the critical accretion luminosity 
from the central star and R is the radius of the star. As a result, the 
magnetosphere is surrounded by a thick (h - r) plasma torus containing cooler, 
denser plasma with a scale height h c << h (see Fig. 1). A central corona and 
mass loss via a wind are expected due to heating of the surface of the disk by 
acoustic flux, magnetic flaring, and radiation from the neutron star .* 3 '! 4 In 
Sco X-l, the existence of radio lobes 15 is direct evidence for mass ejection 
from the system. 

The plasma in the inner disk is expected to be spatially inhomogeneous as 
a result of three mechanisms that can cause clumping. First, vertical 
eruption of magnetic flux carried inward and amplified by the accretion flow 
creates strong spatial inhomogeneities m the density and magnetic field 
strength. 14,1 ® Second, for near-critical accretion rates, the disk- 
magnetosphere boundary lies only a few stellar radii from the surface of the 
neutron star; under these conditions, the inner disk is expected to be 
thermally unstable to the formation of clumps . 17-20 Third, the relative 
motion of magnetospheric and disk plasma drives the Kelvin-Helmholtz 
instability to large amplitudes near the inner edge of the disk, where the 
disk plasma is partially confined by the pressure of the magnetospheric field 
(ref. 21; Aly, Ghosh, & Lamb, reported in ref. 22). 

Assuming that the cooler, denser plasma dominates the mass flux onto the 
neutron star, the disk ends where BpB^S ~ Pc v r r % }l c • 23 Here B p and B 0 are the 
poloidal and toroidal components of the magnetospheric field, & is the width 


of the 

boundary 

layer, p c and v r are 

the 

density and radial velocity 

of 

the 

plasma 

in the 

boundary layer, and 

the 

beat frequency Qg = Q c -n s 

is 

the 
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difference between the angular velocity Q c of the plasma clumps and the 
rotation frequency fi s of the star (see Fig. 2). We expect f) c k Cgo, the 
Keplerian angular velocity at the radius 23 ro n 0 . 52fi 4 / 7 (2GM) -l/7^-2/7 0 f the 
boundary layer. Here p and M are the m,ngnetic moment and mass of the star and 
M is the mass flux through the boundary layer. The spread in beat frequencies 
within the boundary layer is then &n B a j&0(</3r| o S = 1 . 5( S/r 0 ) tlgg . 

In the boundary layer, plasma clumps drift radially inward as they lose 
angular momentum by interaction with each other, with the surrounding gas, and 
with the stellar magnetic field. Plasma is stripped from the clumps by 
interaction with the magnetospheric field and plasma. Important mechanisms 
are likely to include stripping by surface stresses resulting from the Kelvin- 
Helmholtz instability 25 or by reconnection of magnetic flux in the clump with 
the stellar flux, f;) perhaps in a manner analogous to the reported 26 stripping 
of plasma from comets as they move through regions of alternating field 
direction in the solar wind. Once stripped, plasma is quickly brought into 
corotation with the neutron star and accretes to the stellar surface, 
producing X-rays there. 

Unless the stellar field is axisymmetr ic , aligned with the rotation axes 
of disk and star, and centered in the star, the interaction of a given clump 
with the magnetosphere is greater at some stellar azimuths than at others. As 
a result, the rate at which plasma is stripped from the clump is greater at 
some azimuths than others (see Fig. 3}, We refer to the stellar azimuth(s) at 
which the stripping rate is greater for a given clump as the special 
direction(s) for that clump. The pattern of the special directions and the 
steepness of the variation of the stripping rate with azimuth affect the 
harmonic content of the X-ray intensity waveform produced by the clump. For a 
pattern of special directions with k-fold symmetry, the pattern of mass flux 
to the stellar surface and hence the X-ray waveform repeats with frequency fig 
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= kn B (k = 1, 2, 3, ...), Suppose, for example, that the rotation axis of 
the star is aligned with that of the disk but its magnetic axis is tilted with 
respect to its rotation axis For a boundary layer magnetic field that is 
perfectly dipolar, the magnetic pressure on a clump is « l+3cos 2 (Q B t+0) , where 
0 is the phase of the clump. Thus, if the rate at which plasma is stripped 
from clumps is determined entirely by the magnetic stress acting on the 
clumps, k = 2. For the same geometry, reconnection of the magnetospheric 
field previously entrained by the clumps with the local magnetospheric field 
varies with frequency 2fl B , also giving k = 2. However, if there is the 
slightest asymmetry between the two magnetic poles, k = 1, even though most of 
the power in the clump waveform may be at 2fl B . Reconnection to the local 
magnetospheric field of a component of the magnetic field of a clump that 
remains fixed in the corotating frame (due. for example, to the magnetospheric 
interaction with the clump) varies with frequency 0 B . giving k = 1. As long 
as the field in the boundary layer is nearly dipolar, the power in the 
principal peak of the power spectrum typically dominates that at the first 
overtone, due to the smoothness of the dipole field pattern. Even if the 
stripping rate varies as steeply as cos 6 [ (Q B t+0) / 2 ] , the power in the first 
overtone is only 16% of that in the fundamental. 

If the pattern of radiation from the neutron star is not axisyminetric in 
the disk plane at the inner disk, scattering of radiation into the line of 
sight by clumps in the inner disk will create a periodic anisotropy in the 
radiation from the system, causing some modulation of the X-ray intensity seen 
by a distant observer. However, we expect this effect to be much less 
important than modulation of the mass flux to the stellar surface for several 
reasons, including the small solid angle subtended by the clumps at the 
neutron star, the very weak beaming expected from these neutron stars, and the 


6 



suppression of modulation due 


to anisotropic radiation by the plasma 


surrounding the neutron star. The latter two effects are discussed in more 
detail below. 

Power Spectruw 

The theory of random processes (see, for example. Rice 27 ) can be used to 
calculate the power density spectrum of the X-ray intensity time series 1 ( t ) 
produced by the collection of clumps orbiting within the boundary layer. The 
X-ray intensity waveform produced by exactly N clumps may be written 

N 

l N (t) = I Q + E F(t-t j ) , (1) 

1 = 1 

where Iq is a possible constant intensity, F(t;0) is the waveform contributed 
by a single clump, and 0j is the azimuthal phase of the ith clump at the time 
it enters the boundary layer. Quite generally, an arbitrary clump wavefoim 
may be expanded in harmonics of the beat frequency: 

» 

F(t;0) * [ A + Z B n cos (n ( 2nf B t - 0) - a n )] G(t), (2) 

n = l 

where fg - Hg, 2n , a n is the phase of the nth harmonic in the clump waveform, 
and G ( t ) is an envelope function that describes the overall shape of the 
waveform. A, B n , a n , and 0 may he treated as random variables. Equation (2) 
holds even if the special directions are different for each clump, as long as 
they have the same symmetry and remain fixed in the frame of the sta’-. 

The power spectrum produced by the clumps depends on whether they enter 
the boundary layer at random times, whether their positions are correlated, 
and whether stripping of a clump occurs independently of the presence of other 
clumps nearby. For simplicity, we assume initially that correlations are 
absent; later we indicate the effects that such correlations would have. 

Assuming that the number of clumps entering the boundary layer in the 
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observing interval is Poisson distributed, that 0{ is uniformly distributed, 
and that t^ and 0^ are uncorrelated, the expected power density spectrum of 
the time series (1) 1 b 

V 

33 

<P(f)> = 2<l> 2 8( f ) + 2<I> + 2<I> 2 7 RN 2 g(f) + <i> 2 S Y n 2 [g ( f-nf B ) + g(f+nf B )] 

n°l 

(3) 

where <r> = I 0 +<N><A> is the mean intensity, y rn 2 => <NxA 2 >/2<!> 2 , and ? n 2 = 

<N><B n 2 >/4<I> 2 . The structure function g(f) is proportional to |s(f}! 2 , where 
s(f) is the Fourier transform of Gft), but is normalized to unit area. It 
includes the ’lifetime* broadening (FWHM Af L = l/irtl produced by the finite 
duration r of the waveform contributed by a single clump. 

The first term in equation (3) is the power at zero frequency due to the 
mean intensity while the second is the power density in the Poisson noise 
produced by photon statistics. The term proportional to g(f) describes the 
red noise caused by time variations in the photon arrival rate. The presence 
of this additional noise follows immediately from the fact that the mean X-ray 
intensity contributed by each clump is positive and persists only for a finite 
time. All clumps that contribute to the oscillations necessarily contribute 
to the red noise. (The converse need not be true.) The terms proportional to 
g(f-nfg) describe the power density near nf B produced by the oscillating X-ray 
intensity. For most waveforms, the terms proportional to g(f+nfg) may be 
neglected . 

The repetition frequency of the oscillations in the X-ray wave train 
contributed by a single ciump varies in time as the clump drifts inward across 
the boundary layer since its angular velocity varies with radius. This 
variation may be described by introducing a distribution of clump angular 
velocities, In addition, one expects a distribution of clump lifetimes rather 
than a single lifetime. These effects may be included by averaging equation 
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(3) over the expected distribution p(fg.T), which yields equation (31 with g 
everywhere replaced by g * </>g>. In the case of the oscillation terms, g 
equals g n , the convolution of g with />( fg<'n ; t) /n. Thus, if g and p may both 
be approximated by a Lorentzian (or Gaussian) profile, g is a Lorentzian (or 
Gaussian), After this averaging, <P(f)> has peaks at f n - nf B , where f B is 
the frequency </>f B >, with FWHMs that are. loosely speaking, the sum of &fg and 
the width n&fg s ndllg/ 2ir caused by velocity shear in the boundary layer. 
The power in the zero-frequency, red noise, and oscillating components is 


C 0 + f 2 

P 0 * ] <P( f )>df * <I> 2 . 

Jo 


(4) 


Pr N = 2 <I> 27 RN 2 ) e(f) df = <1>e W 

P n - <I> 2 7 n 2 ( g n tf)df as <i> 2 7 n 2 ! B n (f)0f - <I> 2 7 n 2 
Jo J -23 


r : ° 


.2- 2 


and 


(5) 

( 6 ) 


Thus, the rms fractional intensity variation caused by the red noise is Trn = 
!P rn /P 0 ) 3 ' 2 while that due to the oscillation at f n is 7 n = (Pn /p O> 1/2 - For 
comparison, the fractional modulation in the oscillation at f n (defined as 
half the peak-to-peak intensity variation of a sinusoid of equivalent power, 
divided by <!>) is m n = \/ 27 n . 

If F ( t ; 0 ) is a pure sinusoid (B n /0 for only one value of n. say m) and if 
the stripping of each clump occurs independently, A>B m for each clump, since 
then the X-ray intensity contributed by each clump must be positive at all 
times. This implies 7 RN >7 m - If the assumptions concerning tj and made in 
deriving equation (3) are also satisfied, equations (5) and (6) show that P m 
is at most O.SPftft. If F ( t : 0 ) is a pure sinusoid but <pj is not uniformly 
distributed, the oscillation and red noise Fourier amplitudes interfere, 
introducing an additional phase-dependent term =c Re g* ( f )g( f-mf B ) in equation 
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(3) that changes the shape of the spectrum near mfg/2 but alters Pn/^RN by 
only a small amount. Such interference usually also occurs if 0 j is 
correlated with tj. However, for some correlations between 0j and tj this 
interference is absent if t* is uniformly distributed, but the oscillation 
Fourier amplitudes contributed by different clumps interfere constructively 
unlike the red noise amplitudes, causing P m to exceed P^ by a large factor. 
Alternatively, if F ( t ; 0 ) is a pure sinusoid but is scaled by certain functions 
H ( 0 ) , such as cos 2< l(0) for q>>l , P m can approach Prn- Finally, if F(t;0) x 
cos 2c l({2Tfmfgt-0)/2] , with q>3, P m exceeds Pp^* The power density <P(f)> 
used here is related to the power density Pj defined in ref. 28 by the 
expression Pj = <P(f)>/<!>. (Note that although the second of eqs, (1) of 
ref. 28 displays the angular frequency o), Pj in this equation is normalized to 
power per circular frequency interval.) 

Suppose the principal peak in the power spectrum is at the frequency fp ° 
nfg. According to the physical model described above, as the mass flux M 
through the boundary layer changes, the radius rp of the layer changes also, 
causing a change in Ago- Since f s may be treated as constant for times much 
shorter than the timescale for changes in the spin rate of the neutron star, 
here ~ 10 5 years, the change in flgo results in a change in f p . Consider now 
the relation between f p and the X-ray intensity I observed at earth in a given 
energy band. If no mass was lost from the system after passing through the 
boundary layer, if all observed X-rays came from the surface of the neutron 
star, and if there were no changes in the geometry of the emission region or 
the X-ray spectrum, I would be proportional to M, However, tn the bright 
galactic-bulge sources there are strong indications 29-32 that these conditions 
are not satisfied. We therefore consider the more general relation I = 
I E (M/M E )P, where Ip = tiL e in terms of the Eddington luminosity L E = 1.257 
x10 38 (M/Mq) ergs s - ^ for hydrogen plasma. Here q is a conversion factor that 
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taxes into account the distance, as well as geometric and bolometrlc 
corrections. (If I Is the EXOSAT 1 - 1 C keV count-rate. D a Is the distance in 
units of 8 kpc, and if the source emits Isotropically with a spectrum like GX 
5 - 1 . T| » 1 . 42xl0 -3 5p a ‘2 c j ere' 1 ; Cf. ref. 2) Then (cf. ref. 9) 

fp “ n(a(I/I R ) 3/7 P - f s ], (7) 

where a 1 .57xl0 3 (B d /10°Gr 6/ ' 7 (R/10 6 cm)~ 15/7 (M/M Q ) 5/7 . Here B d is the dipole 

component of the stellar magnetic field. A frequently reported quantity 2 - 4 is 

the logarithmic derivative of fp with respect to the X-ray count-rate which. 

in this model, is given by (cf. ref. 9) 

dfln T p 3 1 

a ^ s — , (8) 

dfin I 7(3 (l - w 8 ) 

where w s = n s /n«o ls the fastness parameter. 21 
X-Ray Pulsations 

X-ray pulsations produced by rotation of the neutron stars in the bright 
galactic-bulge X-ray sources have not yet been detected. The canonical view 
is that the magnetic fields of these neutron stars are so weak (B d < 10 R G) 

that they do not channel the accretion flow (see. e.g., ref. 33). In 
contrast, the model of quasiper iodic oscillations developed here requires 
neutron star magnetic fields -10 9 G in the bright bulge sources exhibiting 
such oscillations. Such fields could channel the accretion flow to some 
extent 34 and hence could produce some X-ray beaming. However, in the context 
of the present model there are three distinct physical effects that, for 
distant observers, suppress pulsations at the rotation frequency of the 
neutron star. 

First, X-ray beaming from the neutron star surface is much less in the 
present model than in canonical accretion-powered pulsars. The relatively 
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weak magnetic fluid of the neutron star is less effective In channeling the 
accretion flow.' 1,1 Radiation pressure forces within the magni ..osphcre product^ 
by the high luminosity of the star support the accreting plasma," 15 causing it 
to settle over u large fraction of the star's surfnce. Evidence that 
nccreting plasma falls over a substantial fraction of the surface comes from 
analyses of the hard component In the spectrum, which yield emitting areas 
comparable to the surface area of a neutron star. 4 - 31 ' 52,30 Such a broad 
distribution of accreting plasma and the resulting broad X-ray beam produces 
modulation that is weak and is at relatively high harmonics of the rotation 
frequency of the neutron star, The modulation may be further weakened by 
gravitational bending of the X-ray photon ray paths in the strong 
gravitational field of the neutron star (K. Wood, private communication; see 

ref. 37) . 

Second, the torus and corona around the neutron star (see Fig. 1) tend to 
destroy X-ray pulsations caused by beaming while leaving unaffected the 
quasiperiodic oscillations caused by modulation of the accretion rate. The 
existence of central coronae with dimensions -lo’ 7 cm and optical depths -4-12 
is strongly indicated by analyses of the X-ray spectra of the bright galactic- 
bulge sources . 4 * 32,38 After passing through a corona of electron scattering 
optical depth r es >4 , the observed fractional modulation due to beaming, m 0 |} g , 
is much less than m s tar' wliere m star tlle fractional modulation that would 
be observed in the absence of the corona (for calculations that illustrate 
this behavior in a different context, see Figs. 3 and 5 of ref. 39). In the 
present model, the quasiperiodic oscillations are due to variations in the 
mass flux onto the neutron star, Thus, even if the emission from the neutron 
star were perfectly isotropic, these oscillations would still be observable. 
The fractional modjlation produced by variation of the mass flux is unaffected 
if the mean time for photons to propagate through the corona is less than the 
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oscillation period, i.e., if r cor r es /c < 1 / - p (see Fig. 4). Taking r cor •= 
10 7 cm, r cs = 10, and m star “0.1 as an example, the fractional modulation at 
the rotation frequency would be <<1 whereas the fractional modulation in 30 Hz 
quasiperiodic oscillations would be unaffected, i.e., equal to the fractional 
modulation in the accretion rate. 

Third, the thick torus and disk corona surrounding the neutron star 
prevent us from seeing X-rays coming directly from the neutron star except 
when our 1 ine-of-s ight is close to the star's rufation axis, in which case 
modulation due to beaming is further weakened. The inner disk is expected to 
be both geometrically thick and optically thick to electron scattering (see 
discussion t Fig. 1). There is strong evidence that in bright galactic-bulge 
sources the outer disk is also both geometrically and optically thick. 41- 43 
As a result, we do not see X-rays coming directly from the neutron star except 
in systems with small Inclinations. 33 Assuming that the rotation axis of the 
neutron star is perpendicular to the plane of the disk, such systems are 
viewed from neai the rotation axis, a direction in which modulation due to 
beaming is absent, 

If our analysis of the weakness of X-ray pulsations at the rotation 
frequency of the neutron star is correct, such pulsations should be most 
apparent in sources with relatively low intrinsic luminosities and optically 
thin cororiae. 

Discussion 

The variation of oscillation frequency with X-ray intensity depends on 
the properties of the neutron star and the gross properties of the accretion 
flow, whereas the shape of the power density spectrum is determined by the 
physics of the boundary layer. Consider first the change in oscillation 
frequency . 
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Equation (7) fits the variation in f p with I observed 2 in GX 5-1 for a 
substantial range of neutron star magnetic fields and stellar spin rates, as 
illustrated by Table 1 (see also 2 - 9 ). Models 1-3 assume that f p is the fund- 
amental beat frequency; Models 4-6, that f p is the first overtone. All assume 
that the photons in the so-called ’soft 1 component of the spectrum come from 
the inner disk (see Fig. 1) whereas the photons in the 'hard' component come 
initially from near the surface of the neutron star 12 • 31 - 92 - 36 (but see 4 ). In 
Models 1 and 4 the 40* observed change * 1-10 keV X-ray intensity is assumed 
to be due mostly to a change in the fraction of the radiation from the neutron 
star that is scattered into our line of sight, with the change in the mass 
flux through the boundary layer (and hence the accretion luminosity) only 10*. 
In Models 2 and 5 the change in intensity is assumed to be due mostly to vari- 
ation of the mass flux, which changes by 40*, whereas in Models 3 and 6 the 
change in the mass flux (80*) is much larger than the change in the X-ray 
intensity, as could be the case if most of the photons from the neutron star 
do not reach us. For each value of P, M^ was fixed by assuming that a l-18keV 
EXOSAT count rate of 4850 ct s _1 corresponds to the critical luminosity and 
noting that the count rate 1^ corresponding to is 2427 ct s -1 (we note, 
however, that the value of T| quoted above yields a luminosity at 4850 ct s _1 
of 1 . 9Lg for D = 8 kpc and M = 1.4 Mq). M 2 is then determined by I 2 , which is 
340 ^ ct s _1 . Given the range in mass accretion rate, the dipole component of 
the neutron star magnetic field B<j (=2p/R 3 ) is determined entirely by the 
variation A p of f p from ^ to I 2 . Once is fixed, the stellar rotation fre- 
quency f s is determined by f p . Both B<j and f s are quite sensitive to the 
assumed variation in M. The fastness parameter w s , which is independent of n, 
ranges from 0,95 for Models 1 and 4 to 0.51 for Models 3 and 6. 

Equation (7) can also be fit to the f p -I correlation observed in Cyg X-2 
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as well as that observed in Sco X-l during local minima of the flaring state, 
for similar neutron star magnetic fields and spin rates (see also 4 ' 6 ). The 
currently allowed variation in the parameters among different models of the 
same source is greater than the variation from source to source, for these 
three systems. 

Indirect support for the beat-f requency model is provided by the fact 
that the neutron star rotation rates and magnetic field strengths inferred 
from it are consistent with the idea 44-46 that systems like GX 5-1 and Cyg X-2 
are progenitors of the millisecond rotation-powered pulsars .2 ■ 8 • ® Obviously, 
detection of weak pulsations at the predicted neutron star spin rate would 
provide the most direct evidence in favor of this model. As Table 1 shows, 
measurement of f s would determine the harmonic number n and tightly constrain 
the change in mass flux through the boundary layer. However, it may be 
possible to constrain the beat-frequency model even without direct knowledge 
of f s by extending the fp _ I relation to higher intensity and increasing the 
precision of the measurements, since the shape of the theoretical fp-I curve 
is different for different values of the parameters. 

Consider now the evidence provided by the shape of the power spectrum. 
In the following discussion we assume that the average power spectra reported 
for GX 5-1 and Cyg X-2 accurately reflect the shapes of the individual spectra 
that have been averaged. If this is not the case, the model parameters 
appropriate to these sources will of course be different from those inferred 
here . 

The width of the peak due to quasiperiodic oscillations depends on Af L 
and nAfg whereas the shape of the red noise depends only on Af^. Thus, fits 
of model spectra to observed spectra provide estimates of both quantities. 
Figure 5 compares one spectrum observed in GX 5-1 with two model spectra. In 
GX 5-1, the power in the oscillations is comparable to the p ower in the red 
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noise clown to the l owes t frequencies so far e xplored . when the oscillations 
are most prominent. Recall that if the X-ray intensity waveform produced by 
each clump is sinusoidal and the assumptions leading to equation (3) are 
satisfied, P m can be at most 0.5 Pr n , In the sinusoidal model Pj does indeed 
equal 0.5 Prn, but half the noise power is at frequencies below the lowest 
frequency so far explored, due to the distribution of clump lifetimes. In 
the other model, the waveform is not a pure sinusoid, with the result that P-j 
is I.IPrn. As Figure 5 illustrates, these two models can be adjusted to fit 
the GX 5-1 observations over the frequency range so far studied. For the non- 
sinusoidal model, an acceptable fit is not possible for Af s > 0.5 Hz. The 
fact that the lifetime broadening inferred from the red noise spectrum 
accounts for a substantial fraction of the width of the peak at fp in the 
avera ged power spectra suggests that the widths of the peaks in the i ndividua 1 
spectra are not very different. 

Note that if the clumps are the dominant source of accreting plasma ( I 0 
<< <N><A> ) and <A 2 > - <A> 2 , one can estimate the number of clumps in the 
boundary layer, since then <N> - 1/2 Tr^ 2 . The result is <N a > = <N’b > ~ 30 
for the sinusoidal model and <N> - 120 for the other. Both models can also be 
adjusted to give good agreement with the spectra reported 4 in Cyg x-2. 
Because the spectra they predict have different low-frequency behavior and 
different harmonic content, future observations can distinguish between them. 

Analysis of the harmonic content of the oscillations can help to 
determine whether the observed peak is the fundamental beat frequency or an 
overtone. If f p represents the fundamental (n=l). power may be expected at 
integer multiples of f p but not at f p /2. If instead f p is the first overtone 
(n=2), some power may be expected at the fundamental beat frequency fp/2. The 
distribution of power among the harmonics of f p also constrains the physics of 
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the process by which plasma is stripped from the clumps: the weaker the 
overtones, the more linear the process. The absence of a prominent feature at 
2fp indicates that the variation in stripping rate cannot be as steep as 
COS 2f l [ { 2irf [jt+0) / 2 ] with q>4. 

The ratio of the power Pp at the principal oscillation frequency to the 
power Prn in the red noise depends on the harmonic content and fractional 
modulation of the X-ray waveform produced by a single c lump . The relatively 
large apparent value of this ratio (-1) in GX 5-1 and Cyg X-2 suggests a 
large fractional modulation. If such a large value is confirmed when the 
spectrum is extended to lower frequencies, it would indicate substantial 
azimuthal variation of the conditions in the boundary layer, favoring magneto- 
spheric geometries in which the axis of the stellar field is tilted by a 
substantial angle or is offset from the center of the star by a substantial 
distance. In both sources Pp is observed to vary. If the clump waveform 
remains unchanged as Pp varies. Prjj is proportional to Pp, as observed in GX 
5-1. 2 More generally, Prn need not be proportional to P p and may even 
increase as Pp decreases, as has been reported in Sco X-l 2 and Cyg X-2. 6 If 
the changes in Pp/PRN are due to variations in the shape of the clump 
waveform, the relative strengths of the harmonics of fp must also change, 
whereas if they are due entirely to changes in the fractional modulation of 
the waveform, the relative strengths must remain unchanged. Thus, precise 
measurements of the power spectrum can show which effect is more important. 

The ratios P p /Po find Prn/ p o are related to the rms fractional modulation 
of the waveform produced by all the c lumps involved . The substantial ( ~ 6 — 1 0?o ) 
fractional modulations observed 2 ' 4,6 favor models like the present one in 
which the modulated X-rays come from the neutron star, since that is where the 
bulk of the gravitational energy becomes available.^- 2 If most of the photons 
from the star reach us without gaining or losing appreciable energy, the 
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oscillations and red noise should be most prominent in the spectral component 
that comes from the neutron star. However, photons from the star will suffer 
an energy shift AE a E after only (m e c 2 /kj}T r ) scatterings, if T r > T e , or 
(m e c 2 /l< 0 T e ) scatterings, if T r < T e . 47 Thus, those photons from the star that 

reach us only after scattering through a plasma of optical depth r ~ 10 (for 

T r a 2 keV and T e - 1-7 keV 4 <32,38^ will have energies characteristic of the 
temperature of the plasma rather than that of the surface of the star. As 
shown above, the fractional modulation of the oscillations in the flux of 

these photons is unaffected for a plasma of this optical depth if its 

dimension is < 10 7 cm. In any case, the X-ray spectral dependence of the 
variation in intensity due to the red noise at frequencies not too far below 
fp should be similar to that of the variation due to the quasiperiodic 
oscillations . 

From the narrowness of the observed peak at fp one can show as follbws 
that the lifetime i of a clump ic long compared to the time for it to fall 
freely across the boundary layer, and hence that its angular velocity n c must 
be Keplerian. The time required for a clump to cross the boundary layer in 
free-fall is 21 Atff - ( l/fi«g) • (S/rg) 1 / 2 , where (1 K0 is the Keplerian angular 
velocity at rg. Now Af p > AfL = l/ift. Further, fp ~ (nQ^g/2ir) • ( l-w s ) . Thus, 

r r 0 1/2 2 /n f 

n K 0 r ( — ) » n KO r - 10-200 . (9) 

Atff 8 ^“s * f p 

Given that fi c is closely Keplerian in the boundary layer. Af p > nAfg = 
(n/2n) • |3n c /ar | 0 8 = (3n/4tr) fij<g( S/rg) . Using this relation, one can show 

8 /r 0 < (2/3) ( l-co s ) (Afp/fp) - lO^-lO -2 . (10) 

Noting that v r ~ S/r, one can also show 
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V v 0 < (n/3)(l-U g ) 2 (^f p /fp) 2 - 10" 2 -10' 4 

An independent estimate of t is provided by the noise component, of the 
spectrum, which must flatten below the critical frequency f c ~ 1/2HT. 

The quasiperiodic oscillations seen in Sco X-l are observed to disappear 
altogether at times. 5 - 6 If they are correctly interpreted in terms of the 
present model, the fact that Af p has been observed to increase substantially 
during a transition from flaring to quiescence (Priedhorsky , private 
communication) suggests that conditions within the boundary layer are 
variable. The broadening could be due to increased shear or a shorter mean 
clump lifetime. Fits of equation (3) to the red noise and the oscillation 
peak can indicate which is more important. If conditions in the inner disk 
change sufficiently, the formation of clumps could be Inhibited. We note also 
that if the density and radius of the corona above the inner disk increase to 
the point that r cor r PS /c >> 1 / f p , the oscillations will be destroyed by 
electron scattering in the corona. In this case, changes in the X-ray 
spectrum characteristic of the development of a denser and more extensive 
corona should be observable. Thus, simultaneous measurements of the X-ray 
spectrum and power spectrum can help to determine the relative importance of 
these processes. 

The present model suggests a natural explanation for the 'dips' occasion- 
ally seen- in GX 5-1 when the X-ray intensity is near its minimum value. 31 In 
our model the radius r 0 of the inner edge of the disk is quite close to the 
centrifugal radius r c . Thus, when the mass accretion rate is near its minimum 
value, tq may briefly equal or exceed r c . When this happens, the mass flux 
from the inner edge of the disk to the neutron star temporarily ceases 21 - 34 
and one observes only the 'soft' component from the disk. 31 
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In closing, we emphasize that the model described here may be relevant to 
other galactic-bulge X-ray sources. For example, the 2 Hz oscillations seen 48 
in long flat-top bursts from the Rapid Burster may be a beat frequency 
phenomenon. If so, our model implies that the Rapid Burster has a surface 

magnetic field - 10 11 gauss, assuming f s - 20 Hz. The absence of a so-called 

'soft' component in the X-ray spectrum of the Rapid Burster (V. Tanaka, 
private communication) may also indicate that this source has a significant 
magnetosphere. If the 10“ 3 Hz oscillations reported 48 in the 8 sec accretion- 
powered pulsar 4U1626-67 also prove to be an example of this phenomenon. 9 it 
would demonstrate that the thermally unstable region of the inner disk is not 
essential to the development of clumping, since such a region is absent in 

this source. Our model then gives an estimate of the surface magnetic field 

of - 10* 3 gauss. This estimate is independent of the similar estimate derived 
from the observed spin-up rate, 50 but consistent with it. Finally, we note 
that the present work, when scaled appropriately, may be relevant to some of 
the quasiperiodic oscillations observed in some cataclysmic variables J > 51 > 52 
We thank G. Basinger, M. van der Klis, W.H G. Lewin, and W. Priedhorsky 
for discussing their observations in advance of publication; F. Cordova, S. 
Hayakawa, E.P.J. van den Heuvel, D. Lamb, J. van Paradijs, J. Swank, N. White, 
and G. Zylstra for helpful discussions; and B. Fortner for help in producing 
the figures. This research was supported in part by NSF under grant PHY 80- 
25605 at Illinois and by NASA under grant NSG 7653 at Illinois and grants 
NAGW-567 and NAG 8-497 at Columbia. 
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TABLE 1 


REPRESENTATIVE MODELS OF GX 5-1* 


MODEL 

n Bj 

^s 

Ps 

P 

Mi 

m 2 

r 01 

r 02 

1 

1 5.5 

384 

2.6 

3.6 

0,82 

0.91 

3.1 

3.0 

2 

1 19.9 

88 

11.4 

1.0 

0.50 

0.70 

7.4 

6.7 

3 

1 31.6 

39 

26.0 

0.58 

0.30 

0,55 

11.1 

9.4 

4 

2 12.3 

191 

5.2 

3.6 

0.82 

0.91 

4.9 

4.7 

5 

2 44,6 

44 

22.7 

1.0 

0,50 

0,70 

11 .7 

10.6 

6 

2 70.8 

19 

51.8 

0.58 

0 30 

0.55 

17.6 

14.9 

•All models assume a 

1.4 Mg neutron star of 

radius 

10® cm 

and are 

fit to f p 

=■ 20.07 H 7 . and A p =» 

16.8 Hz, 

The symbols 

and their units are: 

harmonic 

number 

n, dipole component of the 

stellar magnetic field Bj 

(10 9 G). 

stellar 

spin frequency f s (Hz), rotation period P g 

(ms), mass flux 

at 2427 

ct s" 1 Mj (Mg) , mass 

flux at 

3403 ct 

s" 1 M 2 

( Me ) - 

boundary 

layer radius rgj 


{ 10 6 cm) nt and tq 2 (10 6 cm) at M 2 - Mg is the Eddington luminosity of 
the neutron star. 
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FIGURE CAPTIONS 


Fig. 1 Side view of tin* plusnm flow pattern expected around n weakly magnetic 
neutron star .accreting from a disk at near the critical rate. 12 Shown are the 
neutron star (black circle), magnetosphere (dark shaded circle), dense clumped 
plasma (dark shading), torus (moderate shading), and central corona (light 
shading); 6 is the opening angle of the throat of the toruB. Photons (long 
arrows) enter the torus both from the clumped plasma and from the neutron 
star; photons from the neutron star emerging along the axis of the torus are 
scattered by the throat walls and plasma in the throat and corona. Mass is 
lost from the disk and torus via a wind (short arrows). 

Fig. 2 Magnetospher lc (f) a ) and clump (Q c ) angular velocities in the inner 
disk and boundary layer, showing the meaning of the quantities fig, 60 b. r 0 , 

and &• Insidu r co , plasma is in generalized corotation 2 ' 1 with the star; also 
indicated is the centrifugal radius r c at which plasma in the disk corotates 
with the star. The width of the boundary layer has been greatly exaggerated 
for clarity. 

Fig. 3 Schematic drawing of the narrow boundary layer at the inner edge of 
the disk, from which plasma accretes onto the j urface of the neutron star. 
The spiral depicts the trajectory of a given clump, as seen in a frame 
corotating with the star. In this frame, the clump moves with azimuthal 
velocity r(fl c (r)-0 a ) and radial velocity v r . Plasma is stripped from the 
clump continuously, but the rate of stripping depends on the azimuthal 
position of the clump. The sections of the trajectory where the stripping 
ite is high for this particular clump (the ’special directions’ referred to 
In the text) are depictrd with heavy lines whereas those where it is low are 
depleted with light lines. In this example, the frequency of the observed 
oscillations is equal to the fundamental beat frequency Qg = n c (r)-Q s . If 
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instead the pattern of apodal directions has two-fold symmetry, the transfer 
rate is also high along those sections of the trajectory indicated with heavy 
dashed lines and the frequency of the observed oscillations is the first 
overtone of the fundamental beat frequency. 

Fic. 4 Fractional modulation In quasiperiodlc oscillations seen by a distant 
observer as a function of coronal radius r cor and oscillation frequency f p , 
tor r es « 10. Heavy curve: uniform sphere; light curve: thin shell. For the 
example given in the text, log(r cor fp/c) a -2. (Adapted from ref. 40.) 

Fig. 5 Top: Power spectrum observed in GX 5-1 at 2277-2486 counts s -1 , when 
fp o 20.07 Hz and Af p « 4.2 Hz. 2 The spectrum to the right of the peak has 
been smoothed, A visual estimate of the la error at 2f p is shown. Bottom: 
Power spectrum calculated for a model with F «: cos 6 [(2fifgt+$)/2] and a single 
clump lifetime (solid curve: 7 j = 2,572 3 15 ?3 * 1 • 1 Yrn ” 6.7xl0 -2 , ° 
20.07 Hz, r a 0.085s, Af Si ° 0.47 Hz; thc relative widths Aln(f) of the peaks 
at f lt 2f i , and 3fj are 0.21, 0,12, and 0,085, respectively) and a model with 
F < l+cos(2nf B t+0) and two clump lifetimes (dashed curve: 7 RNa = 7RNb ° \/27i a 
= >i272b = 6 . 3xl0“ 2 , fj - 20.07 llz, T a = 50 s . r b = 0.085s. Af Sa = Af sb 2.6 
Hz). Both models assume the assumptions leading to equation (3) are 
satisfied, a clump waveform G(t) = 0{t)exp(-t/r) . and a horentzian angular 
velocity distribution, which gives gi(f) “ 2r [ [ l+-(2trfTj ) and & n i(f-f n ) 3 
2 t^( l+«nAf S Ti ) • { [2tr ( f-f n )T t ) ] 2 + [ 1+nnAfgTi] 2 )' 1 . Note the marked difference in 
the low-frequency behavior of the two calculated spectra and the much smaller 
difference in harmonic content. 
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